Gravitational Radiation Instability in Hot Young Neutron Stars 
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We show that gravitational radiation drives an instability in hot young rapidly rotating neutron 
stars. This instability occurs primarily in the I = 2 r-mode and will carry away most of the angular 
momentum of a rapidly rotating star by gravitational radiation. On the timescale needed to cool 
a young neutron star to about T = 10 9 K (about one year) this instability can reduce the rotation 
rate of a rapidly rotating star to about 0.076Qk, where Qk is the Keplerian angular velocity where 
mass shedding occurs. In older colder neutron stars this instability is suppressed by viscous effects, 
allowing older stars to be spun up by accretion to larger angular velocities. 



PACS Numbers: 04.40. Dg, 97.60. Jd, 04.30. Db 

Recently Andersson Q| discovered (and Friedman and 
Morsink confirmed more generally) that gravitational 
radiation tends to drive the r-modes of all rotating stars 
unstable. In this paper we examine the timescales asso- 
ciated with this instability in some detail. We show that 
gravitational radiation couples to these modes primar- 
ily through the current multipolcs, rather than the usual 
mass multipoles. We also evaluate the effects of internal 
fluid dissipation which tends to suppress this instabil- 
ity. We find that gravitational radiation is stronger than 
viscosity in these modes and so this instability severely 
limits the rotation rates of hot young neutron stars. We 
show that such stars can spin down by the emission of 
gravitational radiation to about 7.6% of their maximum 
rotation rates on the timescale (about one year) needed 
to cool these stars to 10 9 K. 

The r-modes of rotating barotropic Newtonian stars 
are solutions of the perturbed fluid equations having (Eu- 
lcrian) velocity perturbations 



Sv = am (— 
\R 



(1) 



where R and f2 are the radius and angular velocity of the 
unperturbed star, a is an arbitrary constant, and Y\ m is 
the magnetic type vector spherical harmonic defined by 



?& = [ia + l)]- 1/2 rVx (rVY lm ) 



(2) 



Papaloizou and Pringle || first showed that the Euler 
equation for r— modes determines the frequencies as 
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Further use of the Euler equation (as first noted by 
Provost, Berthomieu and Rocca Q) in the barotropic 
case (a good approximation for neutron stars) determines 
that only the I = m r-modes exist, and that Sv must have 
the radial dependence given in Eq. (0). These expressions 



for the velocity perturbation and frequency are only the 
lowest order terms in expansions for these quantities in 
powers of f2. The exact expressions contain additional 
terms of order Q 3 . 

The lowest order expressions for the (Eulerian) density 
perturbation Sp can also be deduced from the perturbed 
fluid equations (Ipser and Lindblom S): 
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where S^(r) is proportional to the gravitational potential 
(5$ and satisfies 
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Eq. (0) is the complete expression for Sp to order f2 2 . 
The next order terms are proportional to Q . 

Our interest here is to study the evolution of these 
modes due to the dissipative influences of viscosity and 
gravitational radiation. For this purpose it is useful to 
consider the effects of radiation on the evolution of the 
energy of the mode (as measured in the co-rotating frame 
of the equilibrium star) E: 



E= l - 
2 
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This energy evolves on the secular timescale of the dis- 
sipative processes. The general expression for the time 
derivative of E for a mode with time dependence e wt and 
azimuthal angular dependence e %mv is 
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The thermodynamic functions r\ and C that appear in 
Eq. (0) are the shear and bulk viscosities of the fluid. The 
viscous forces are driven by the shear Sa a b and expansion 
5a of the perturbation, defined by the usual expressions 



Scr a b = k{^a5v b + V b 5v a - §6 ab V c 6i> c ), 



5a = \7 a 6v a . 



(8) 
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Gravitational radiation couples to the evolution of the 
mode through the mass 6Di m and current 8 Ji m multipole 
moments of the perturbed fluid, 



6D lm = / dpr l Y; m d 3 x, 



(10) 
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r^pSv + Sp^-Y^x, (11) 



c\l l + l 
with coupling constant 
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The terms in the expression for dE/dt due to viscosity 
and the gravitational radiation generated by the mass 
multipoles are well known ||. The terms involving the 
current multipole moments have been deduced from the 
general expressions given by Thorne Q . 

We can now use Eq. (]?]) to evaluate the stability of the 
r-modes. Viscosity always tends to decrease the energy 
E, while gravitational radiation may either increase or 
decrease it. The sum that appears in Eq. ([?]) is positive 
definite; thus the effect of gravitational radiation is de- 
termined by the sign of w(w + mfi). For r-modes this 
quantity is negative definite: 
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Therefore gravitational radiation tends to increase the 
energy of these modes. For small angular velocities the 
energy E is positive definite: the positive term \5v\ 2 in 
Eq. (|J) (proportional to ft 2 ) dominates the indefinite 
term (Sp/p — 5Q)5p* (proportional to ft ). Thus, gravi- 
tational radiation tends to make every r-mode unstable 
in slowly rotating stars. This confirms the discovery of 
Andersson Q and the more general arguments of Fried- 
man and Morsink Q . To determine whether these modes 
are actually stable or unstable in rotating neutron stars, 
therefore, we must evaluate the magnitudes of all the dis- 
sipative terms in Eq. (^) and determine which dominates. 



Here we estimate the relative importance of these dis- 
sipative effects in the small angular velocity limit using 
the lowest order expressions for the r— mode Sv and Sp 
given in Eqs. ([!]) and (^). The lowest order expression 
for the energy of the mode E is 



E = 
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The lowest order contribution to the gravitational radi- 
ation terms in the energy dissipation comes entirely from 
the current multipole moment SJu. This term can be 
evaluated to lowest order in ft using Eqs. (Ill) and (ITT 
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The other contributions from gravitational radiation to 
the dissipation rate are all higher order in Q. The mass 
multipole moment contributions are higher order because 
a) the density perturbation Sp from Eq. (^) is propor- 
tional to fl 2 while the velocity perturbation Sv is propor- 
tional to f2; and &) the density perturbation Sp generates 
gravitational radiation at order 21 + 4 in w while Sv gen- 
erates radiation at order 21 + 2. 

The contribution of gravitational radiation to the 
imaginary part of the frequency of the mode 1/t gr can 
be computed as follows, 



T~GR 




(16) 



Using Eqs. (0)-(0) we obtain an explicit expression for 
the gravitational radiation timescale associated with the 
r-modes: 
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The time derivative of the energy due to viscous dissi- 
pation is driven by the shear Sa a b and the expansion 5a 
of the velocity perturbation. The shear can be evaluated 
using Eqs. (Q) and (||) and its integral over the constant r 
two-spheres performed in a straightforward calculation. 
Using the formulae for the viscous dissipation rate Eq. (Q) 
and the energy Eq. (|L4|), we obtain the contribution of 
shear viscosity to the imaginary part of the frequency of 
the mode, 
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The expansion 5a, which drives the bulk viscosity dis- 
sipation in the fluid, can be re-expressed in terms of the 



2 



density perturbation. The perturbed mass conservation 
law gives the relationship da — —i{uj + mQ)Ap/ p, where 
Ap is the Lagrangian perturbation in the density. The 
perturbation analysis used here is not of sufficiently high 
order (in fl) to evaluate the lowest order contribution to 
Ap. However, we are able to evaluate the Eulerian per- 
turbation 5p as given in Eq. (Q). We expect that the 
integral of \5p/p\ 2 over the interior of the star will be 
similar to (i.e., within about a factor of two of) the inte- 
gral of \Ap/p\ 2 . Thus, we estimate the magnitude of the 
bulk viscosity contribution to the dissipation by 

j, + / C W,3^ (ig) 

t b 2E J p 2 

Using Eqs. (f|) and (|l|) for Sp/p and E, Eq. (|l|) becomes 
an explicit formula for the contribution to the imaginary 
part of the frequency due to bulk viscosity. 

To evaluate the dissipative timescales associated with 
the r-modes using the formulae in Eqs. (|T^)-(|l9|), we 
need models for the structures of neutron stars as well 
as expressions for the viscosities of neutron star matter. 
We have evaluated these timescales for 1.4M Q neutron 
star models based on several realistic equations of state 
. We use the standard formulae for the shear and bulk 
viscosities of hot neutron star matter |9| 

j] = 347p 9/4 T~ 2 , (20) 

C = 6.0x 10- 59 p 2 (uj + mn)- 2 T (i , (21) 

where all quantities are given in cgs units. The timescales 
for the more realistic equations of state are comparable 
to those based on a simple polytropic model p = np 2 
with k chosen so that the radius of a 1.4M star is 12.53 
km. The dissipation timescales for this polytropic model 
(which can be evaluated analytically) are f GR = —3.26s, 
f v = 2.52 x 10 8 s and f B = 6.99 x 10 8 s for the fiducial val- 
ues of the angular velocity fl — y/irGp and temperature 
T = 10 9 K in the I = 2 r-mode. The gravitational radi- 
ation timescales increase by about one order of magni- 
tude for each incremental increase in while the viscous 
timescales decrease by about 20%. 

The evolution of an r-mode due to the dissipative ef- 
fects of viscosity and gravitational radiation reaction is 
determined by the imaginary part of the frequency of the 
mode, 

1 J_ 1 /10 9 K\ 2 




Eq. ( |22] ) is displayed in a form that makes explicit the an- 
gular velocity and temperature dependences of the var- 
ious terms. Dissipative effects cause the mode to decay 



exponentially as e _t / T (i.e., the mode is stable) as long 
as r > 0. From Eqs. Jl7|)-(|l9|) we see that f v > and 
t b > while f GR < 0. Thus gravitational radiation drives 
these modes towards instability while viscosity tries to 
stabilize them. For small f2 the gravitational radiation 
contribution to the imaginary part of the frequency is 
very small since it is proportional to £l 2l+2 . Thus for 
sufficiently small angular velocities, viscosity dominates 
and the mode is stable. For sufficiently large fl, however, 
gravitational radiation will dominate and drive the mode 
unstable. It is convenient to define a critical angular ve- 
locity fl c where the sign of the imaginary part of the fre- 
quency changes from positive to negative: 1/t(£7 c ) = 0. 
If the angular velocity of the star exceeds il c then grav- 
itational radiation reaction dominates viscosity and the 
mode is unstable. 
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FIG. 1. Critical angular velocities for a 1.4M© polytropic 
neutron star with (solid) and without (dashed) bulk viscos- 
ity. Also the evolution of a rapidly rotating neutron star 
(dash-dot) as the star cools and emits gravitational radiation. 

For a given temperature and mode I the equation for 
the critical angular velocity, = 1/t(£! c ), is a polyno- 
mial of order I + 1 in £l 2 , and thus each mode has its own 
critical angular velocity. However, only the smallest of 
these (always the I = 2 r-mode here) represents the crit- 
ical angular velocity of the star. Fig. [l] depicts the criti- 
cal angular velocity for a range of temperatures relevant 
for neutron stars. The solid curve in Fig. |l| represents 
the critical angular velocity for the polytropic model dis- 
cussed above. Fig. || depicts the critical angular velocities 
for 1.4M Q neutron star models computed from a vari- 
ety of realistic equations of state ||. Fig. || illustrates 
that the minimum critical angular velocity (in units of 
y/nGp) is extremely insensitive to the equation of state. 
The minima of these curves occur at T w 2 x 10 9 K, with 
£! c w 0.043 v / 7rGp. The maximum angular velocity for 
any star occurs when the material at the surface effec- 
tively orbits the star. This 'Keplerian' angular velocity 
Hk is very nearly ^^/irGp for any equation of state. Thus 
the minimum critical angular velocity due to instability 
of the r-modes is about 0.065r2x for any equation of state 
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FIG. 2. Critical angular velocities of realistic 1.4M© neu- 
tron star models. 
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To determine how rapidly a young neutron star is al- 
lowed to spin after cooling, we must compare the rate it 
cools with the rate it loses angular momentum by grav- 
itational radiation. We approximate the cooling with a 
simple model based on the emission of neutrinos through 
the modified URCA process p] ]. We compute the time 
evolution of the angular velocity of the star by setting 
dJ/dt = J/t, where J is the angular momentum of the 
star and r is the timescale given in Eq. (^). The result 
is a simple first order differential equation for £l(t) which 
we solve for initial angular velocity f2 = Qk and initial 
temperature 10 11 K. The solution is shown as the dash- 
dot line in Fig. |l|. The gravitational radiation timescale is 
so short that the star radiates away its angular momen- 
tum almost as quickly as it cools. The angular velocity 
of the star decreases from to 0.076f2if in a period 
of about one year [{L2|| . Thus, we conclude that young 
neutron stars will be spun down by the emission of grav- 
itational radiation within their first year to a rotation pe- 
riod of about 13P m in, where P m in = 2it/£Ik- The Crab 
pulsar with present rotation period 33ms and initial pe- 
riod 19ms (based on the measured braking index) rotates 
more slowly than this limit if P m ; n < 1.5ms. 

Our analysis here is based on the assumption that a 
young hot neutron star may be modeled as a simple or- 
dinary fluid. Once the star cools below the superfluid 
transition temperature (about 10 9 K) the analysis pre- 
sented here must be modified jl3| . We expect the r-mode 
instability to be completely suppressed (with Q c = O^-) 
when the star becomes a superfluid fTi)] . This makes it 
possible for old recycled pulsars to be spun up to large 
angular velocities by accretion if they are not re-heated 
much above the Eddington temperature of 10 7 K in the 
process. If non-perfect fluid effects enter above 10 9 K, 
however, the spin down process may be terminated at 
a higher angular velocity than the 0.0760^ figure com- 
puted here. The detection of a young fast pulsar would 
provide evidence for such non-perfect fluid effects at these 
high temperatures. 
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